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Abstract 



The questions of justification of the Gibbs canonical distribution for systems with 
elastic impacts are discussed. A special attention is paid to the description of prob- 
ability measures with densities depending on the system energy. 

Gibbs Distribution 

Let x = (xi, . . . , x n ) be generalized coordinates, y = (yi, . . . , y n ) be conjugate 
canonical momenta of a Hamiltonian system with n degrees of freedom and a sta- 
tionary Hamiltonian H(x, y, A), where A = (Ai, . . . , A m ) are some parameters. Ac- 
cording to Gibbs [1], a distribution with the probability density 



where c = const > 0, (3 = (T is absolute temperature, k is Boltzmann constant) 
plays a key role in statistical consideration of Hamiltonian systems. The constant c 
is chosen due to the normalization condition of density p. 

Given an invariant measure with the density (1), we can introduce an mean 
energy 



and average generalized forces (constraint reactions A = const), corresponding to 
the parameters A: 



p = ce 



-(3H 



(1) 
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1 < i < m . 



(3) 
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Relations Li = fi(/3, A) are considered as equations of state. 
As it was shown by Gibbs, 1-form of heat gain 



dE + J^A+dK ( 4 ) 



satisfies the axioms of thermodynamics: the form (3uj is exact (/3u = dS, where S((3, A) 
is the entropy of a thermodynamical system). In particular, the form ui is an exact 1- 
form under fixed values of (3. Thus, according to Gibbs, to any Hamiltonian system 
(provided that the integrals (2) and (3) exist and depend smoothly on A and 0) there 
can be associated a thermodynamic system with external parameters Ai, . . . , A m , 
the internal energy (2), and the equations of state (3). The relations (2) and (3) 
can be simplified by introducing statistical integral 

Z(0, A) = J e^ H d n xd n y. (5) 

Hence, 

dlnZ IdlnZ 

and therefore (3u = dS, where 
Thermodynamics of Billiards 

Billiard is a mass particle performing the inertial motion in domain D of three- 
dimensional Euclidean space and reflecting elastically from its boundary dD. We 
can consider a more general case when there are n-identical particles in the domain D 
not interacting with each other (in particular, not colliding with each other). Such 
a system is a universally recognized model of rarefied perfect gas. 

Let qi = (x{, yi, Zi) be a set of Cartesian coordinates of the i-th particle of unit 
mass with momentum pi = (ij, y*, i«). Dynamics of the system in the domain D is 
defined by a Hamiltonian 

^ 2 

Since this function does not contain any information about the geometry of the 
domain D, equations (2) and (3) are not applicable immediately. In this case, one 
can apply the following procedure: statistical integral (5) is calculated first, and 
then relations (6) are used. In our case 



/ / 



3n 
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where v is the volume of D. Therefore, the only external parameter A is the volume v; 
a conjugate variable A is the gas pressure p inside D. Taking into account (8), 
from (6) we obtain known equations of a perfect gas 

^ 3 lrT1 nkT . . 

E=-kT, p= . (9) 

2 v 

Billiards, being systems with one-way constraints, are idealization of ordinary 
mechanical systems with smooth Hamiltonians. When a particle hits the wall, 
the wall deforms giving rise to great elastic forces which push the particle back 
into D. These elastic forces are modeled by potential V u (q). It equals zero in D 
and — outside D. Here / is a smooth function that defines the boundary equa- 
tion dD: f(q) = 0. The large constant v plays a role of elasticity coefficient. It is 
assumed that the boundary does not contain critical points of the function /; in par- 
ticular, boundary 3D is a smooth regular surface. As was shown in [2], as v — > oo, 
solutions of a system with the Hamiltonian 

H =j + V M ( 10 ) 

tend to the motions of a system with elastic reflections in D. 

Application of the Hamiltonian (10) gives corrections to the expression of statis- 
tical integral which depend on the area a of the boundary of D. Thus, the area a 
should be as an external parameter of the perfect gas as a thermodynamic system; 
pressure will be the function of not only volume and temperature, but also of the 
surface area of a vessel. 

The meaning of the correction is that the volume v in (8) is replaced by 



provided that / does not have critical points outside D. Taking this fact into ac- 
count, the equation of internal energy E remains the same and the state equation (9) 
is replaced by 

nkT ixk , , 

P= ~ r-r ' X= V- (12) 

Since a is a new thermodinamical parameter, we should introduce a conjugate 
variable 

1 din Z _ nkT^fxT 
V ~P^~ v + VxTa' { } 

The relations (12) and (13) constitute a total system of the state equations. 
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Let us indicate the deduction of the formula (11). To do so, we use an obvious 
formula 

/ e~ pVv d 3 q = v+ / e~^d 3 q . 
Jr 3 Jf>0 

According to the saddle-point method, the basic contribution to the asymptotics of 
the second integral as v — > oo is made by the critical points of the potential V. In 
accordance with the assumption, df ^ for / > 0. Consiquently, a set of critical 
points coincides with the boundary dD = {/ = 0}. A non-isolation of the criti- 
cal points results in a certain difficulty under usual application of the saddle-point 
method. Let us pass (locally) into a neighdourhood of the boundary to semigeodesi- 
cal coordinates U\, u 2 , u 3 , where / = Ui [3]. In these variables, the Euclidean metric 
is written in the form 

du\ + a du\ + 2b du 2 du 3 + c du\ , 

where a, b, c are smooth functions of u s . In these variables the desired integral is 
replaced asymptotically by the integral 

I giuije 2 du ± , (14) 

Ju!>0 



where 



VG du>2 du 3 , G = ac — b 2 > . 



dD 



Then, with the help of a standard method [4] , we obtain the asymptotics of the 
integral (14): 



Note now that g(0) = a. 



Probability Distribution 

Now a rigorous deduction of the Gibbs distribution is given only for the case of 
vanishing interaction of individual subsystems. A classical Darwin-Fauler approach 
represents an asymptotical (as n — > oo) deduction of Gibbs distribution from the 
general principles of dynamics in the assumption of the ergodic hypothesis. As it 
is observed by A. Ya. Hinchin [5], this approach repeats in fact the previous mathe- 
matical results, connected with the limiting probability theorems. 

In [6] there suggested another deduction of distribution (1). It is based on the 
fact that the probability density is a single- valued first integral [1]. With the help 
of Poincare method, the conditions, under which motion equations of interacting 
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subsystems do not admit integrals of C 2 -class, independent of energy integral, are 
indicated. These conditions are constructive and, obviously, less strong than the 
assumption of ergodicity. Moreover, a natural Gibbs postulate about thermody- 
namical equilibrum of subsystems under vanishing interaction is used in [6]. 

A statistical analogue of this argument is the deduction of a normal distribution, 
suggested by Gauss. He does not use the central limiting theorem, but the postulate 
that a sample mean is an estimate of maximum of probability at the finite number 
of observations n > 3 (see [7], [8]). 

In connection with the above-said it is usefull to set in order the hierarchy of 
Hamiltonian dynamical systems with respect to the degree of their arbitrariness. 
Let us fix the phase space P of dimention 2n > 4 with analitical structure and 
introduce into consideration a set 7i of all Hamiltonian systems on P with analytical 
Hamiltonians. Certainly, it is supposed that the property of Hamiltonians to be 
analitical on P is concordant with analytical structure of P itself. 

We introduce a sequence of embedded into each other sets of H: 

M G £ GT G K° G K} G . . . G 1C°° G A. (15) 

Here, M, £ and T are the set of systems, which respectively possess the properties 
of intermixing, ergodicity and trasitivity on the energy (2n — l)-dimentional sur- 
faces. Further, JC S is the set of systems, which do not admit the first integrals of 
smoothness class of C S (P) not depending on the energy integral. In addition, the 
case s = corresponds to the continuous integrals: they are locally unstable on the 
surfaces of the level of energy integral and take equal values on the trajectories of 
the Hamiltonian system. The symbol A denotes Hamiltonian systems, which do not 
admit an additional analytical integral. 

One can deduce an analogous chain of embedded sets for the systems with elastic 
reflection as well. 

First of all we should make sure that the neighbouring sets in the chain (15) do 
not coincide with each other. The inequalities Ai^£,£^T,T^JC° can be much 
easily demonstrated by the examples of area preserving mappings of T of the two- 
dimentional torus T 2 = {x, y mod 2n}. Such mappings can be treated as Poincare 
mappings of the energy manifolds cuts of Hamiltonian systems with two degrees of 
freedom. A classical example of mixing transformation drives an automorphism of 
a torus, given by a uni-modular matrix 

"2 l" 
1 1 

The shifts x — > x + a, y — > y + b, where numbers a, b and 2ir are rationally in- 
commensurable, provide us with known examples of ergodic, but not mixing, trans- 
formations. Thus, M. 7^ £. It is considerably more difficult to give examples of 
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transitive, but not ergodic, transformations with an invariant measure. For the 
first instances of such transformations we cite the work by L. G. Shnirelman (1930) 
and A. Bezikovich (1937). They considered continuous automorphisms of a circle. 
Smooth modifications of such transformations are indicated in [9]. 

To proof the inequality T ^ JC° we use an example of transitive area preserv- 
ing transformation T of the square K 2 which leaves the points on its boundary 
immovable. Such an example is built by Oxtoby [10] with the help of theory of 
categories of sets. Let us take four such squares and form one square of quadrupli- 
cated area out of the four (see Figure 1). Identifying opposite sides, we will obtain 
two-dimentional torus, where the mapping T: K 2 — > K 2 is naturally prolonged to 
the contituous area- preserving mapping T: T 2 — > T 2 . We need hardly mention that 
this transformation will no longer be transitive. Still it does not admit non-constant 
continuous integrals. It would be interesting to provide an analytical example of the 
transformation from the set of K,° \ T . 

Inequalities IC k ^ IC k+1 (k = 0, 1, . . . , oo) and /C°° ^ A are derived from the 
results of [11] (see also [12]), where the examples of analytical Hamiltonian systems, 
not possessing additional integrals of C k (C°°) -class, but at same time not admitting 
integrals of C k+1 (C w ) -class, are indicated. 

Let us consider one of the links of the chain of inclusions of (15), say, T C K°. 
The question is, which of the two sets is more massive: T or K,°\T . Apparently, the 
second. However, the answer to this question (as well as its formulation) depends 
on the introduced topology in the space K°. Analogous assumptions are probably 
valid for any pair of the neighbouring sets in (15). 

Classes of systems from (15) may be laid out into a wider class of systems, 
which do not admit additional single-valued complex-analytic first integrals. An 
obstacle to the existence of single-valued holomorphic integrals is the branching of 
the solutions of Hamiltonian systems in the plane of complex time. The discussion 
of this range of questions one can find in the work [13]. 

If we remain within the real examination, then the class A admits a natural 
extension for the dynamics of natural machanical systems. They are decribed by the 
Hamiltonians of the form H = T+V, where T is a kinetic energy, a positively defined 
quadratic form with respect to the momenta, and V is a potential energy, a function 
on the configuration space. All known integrals of such systems are polynomials in 
momenta with single- valued, coefficients on configuration space, (or functions of such 
polynomials). In analytical case, these coefficients are also represented by analytical 
functions. We can show that the existence of an additional polynomial integral of 
the system with the Hamiltonian H = T + V is equivalent to the existence of an 
integral of the system with the Hamiltonian H = T + eV (e is a small parameter) 
as the series in terms of powers of e. 

This problem is more simple and since Poincare times there have been proposed 
efficient methods for its solution [13]. The existence conditions of additional poly- 
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nomial integral of a plane billiard are obtained with the help of complex variable 
function [14]. 

The issue of whether a certain Hamiltonian system belongs to the class A is more 
complex. But essential advances have been made in this field as well, especially for 
the case with small number of degrees of freedom [13]. Difficulties become much 
more severe as we move towards the beginning of the chain (15). Thus, according to 
Kac [15], an efficient verification of the ergodic property of a dynamical system is a 
nearly hopeless problem. Moreover, in many important cases, from the application 
viewpoint, ergodic hypothesis is refuted by the results of KAM theory. For instance, 
as it was established by Lazutkin [16], a billiard inside a plane convex curve (of C 2 - 
class of smoothness) is not ergodic. It does not even possess the transitive property. 
Lack of ergodicity in spatial case was proved in [17] under some additional conditions. 
These examples are directly related with the deduction of Gibbs distribution for the 
perfect gas. 

For small perturbation of an integrable Hamiltonian system with two degrees 
of freedom, Kolmogorov tori cut a three-dimentional energy surfaces. Therefore, a 
perturbed system can no longer be transitive. On the other hand, as it was noted by 
Arnold, such systems admit a nonconstant continuous integral that takes constant 
values in slits between Kolmogorov tori. It's not quite clear yet whether such systems 
have locally nonconstant continuous first integrals which are not identically constant 
in any neighbourhood of every point of the energy surface. A simpler problem 
is whether perturbed systems of general kind with two degrees of freedom admit 
nonconstant integrals of C 1 -smoothness class. 

For systems with n > 3 degrees of freedom, the slits between Kolmogorov tori 
form a connected set everywhere densely filling a five-dimensional energy manifolds. 
Therefore, a principal possibility of the appearance of transitive property arises. This 
is one of the exact statements of the known hypothesis of diffusion in perturbated 
multidimensional Hamiltonian systems. For the purpose of statistical mechanics 
this diffusion hypothesis can be formulated in a less restricted fashion: is it true 
that under great n a perturbed Hamiltonian system of general form does not admit 
nonconstant continuous (or even smooth, of C 1 -class) first integrals on (2n — 1)- 
dimensional energy surfaces? In fact, it is sufficient that this property appeared 
under a small fixed value of perturbing parameter e and a great value of n of weakly 
interacting subsystems. 

Generalized entropy 

Our observations described in previous Section result in a natural assumption that 
the density of probability distribution p is a function of H. The question is: what 
makes Gibbs distribution different from all other distributions of this kind? 

Let z — * f(z) be a nonnegative real function of one variable, /' be its derivative. 
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Following Gibbs, we will consider probability density 

f(PH) 
P jf{(3H)dxdy 



(16) 



assuming that the integral converges over the whole phase space. Here again f3~ l = 
kT. When / = ce~ z , c = const ^ 0, we shall obtain Gibbs distribution. We 
could consider a more general case, when the function / depends also on external 
parameters A (as well as the Hamiltonian H). But we shall not follow this case. 

Let us calculate an average energy E and generalized forces A using (2) and (3), 
with density p determined by (16). Then we can compose 1-form of heat gain in 
accordance with (4). Using direct calculations we can prove 

Theorem. The form uj satisfies axioms of thermodynamics if 

I w/ dxdy I w/ dxdy = Iw/ dxdy I w/ dxdy - (17) 

for all 1 < i, j < m and 

/f OH f OH f 

H fdxdy J — fdxdy = J —fdxdy J H fdxdy (18) 

for all 1 < i < m. 

It is obvious that for the function f(z) = ce~ z these conditions are met. Equali- 
ties (17) and (18) can be rewritten as follows 

OF OF 

A J ^ = A J — , (l<i,3<m), (19) 
F OF OF 

Jaxr- K 'W- (1£8£m) ' (20) 

where 



F = J f((3H)d n xd n y. 



By analogy with Gibbs case, the function F can be called a generalized statistical 
integral. 

From (19) and (20) follows the existence of function x((3, Ai, . . . , A m ), such that 

x OF OF 
A - = "9A-' E = *Jp- (21) 

Therefore, the form of heat gain takes the form 
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Axioms of thermodynamics impose constraints on the form of function x. From (19) 
we obtain a series of inequalities 



8x 8F dxdF n .„ . . 



and the equation (20) yields relations 

8k dF 8x 8F 
3/3 3A~ ~ d\i~dp 



, (1 < i < m) . (23) 



Equalities (22) and (23) denote that functions x and F are dependent. Therefore, 
we can write that x = x(F), at least locally. 

Let <3> be antiderivative of x(-). Then equalities (21) take a simpler form 

Ai = -/*aV £= ar (24) 



Hence 



d ^dp)-d0 d0 -22d^ dXi = 

3/3 



The function 



<9$ 

* = ^-» (25) 

is called an entropy in thermodynamics. 

The form of this function suggests that Legendre transform over (3 should be 
applied. Assuming that 

from the second relation of (24) we will obtain j3 as a function of E and A. We will 
assume E, Ai, . . . , X m independent variables. Then S = S(E,X) and from (25) we 
will obtain potential form of basic thermodynamic relations (24): 
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The Perfect Gas 



Let us apply relations from previous Section to the perfect gas inside domain D of the 
three-dimentional Euclidean space; let v be the volume of D. Remembering that the 
perfect gas is a totality of n equal and not interacting particles performing the inertial 
motion inside D and reflecting elastically from its boundary dD. When taking into 
account arbitrarily small interaction of particles, we will obtain a system without 
additional integrals and therefore we can consider that the density of probability 
distribution is a function of total energy. Let particle interaction tends to zero; 
then we will obtain simple equations for average energy and state equations; these 
equations define thermodynamics of a simplified system, i. e. the perfect gas. Let 
particle mass be equal to unit. Hence, the Hamiltonian for the perfect gas will be 
determined by the following equation 

„2 



2 7 

where Pi = j/j, ij) is momentum of the z-th particle; let qi = (xj, i/i, Zi) be its 
Cartesian coordinates. 

The formula for internal energy has the form 



E = 



V/ d 4£p, 2 /(!Ep?) d 3n pd^ q 



/r*.L»/(§£p?) d 3n Pd 3n 1 
It is independent of volume v: 



where 



b= I f[^> uf)d 3n u 



Variables p and u are connected by simple relations: U; L = y/j3pi. 

Assuming for simplicity 3n = m + 2, we will pass from u\, . . . , u m+ 2 to spherical 
coordinates r, 9i, . . . , 9 m , ip using the following equations 

u± — r cos $i , 

U2 = r sin 9\ cos #2 , 

u 3 = r sin 9\ sin 9 2 cos 9 3 , 



u m = r sin 9\ sin 9 2 . ■ ■ sin 9 m -\ cos 9, 
u m+ i = r sin 9\ sin 9 2 . . . sin 9 m cos tp , 
u m +2 = t sin 9i sin 9 2 ■ ■ ■ sin 9 m sin ip . 
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Here r > 0, < 8j < n (1 < j < m) and ip mod 2tt is an angular coordinate. 
In the new coordinates 



b= '"" +1 /(-J* / ... / / (sm6 1 ) m (sm6 2 ) m - l ...(sme m )d6 1 ...de m d V = 



where T is Euler's gamma-function. By analogy, 



r 



o l+f /-oc m+3 /r 2, 

"" ' " ' -f -TV )<''■■ (28) 



(l + f ) 



2 " V 2 



Now we calculate generalized statistical integral: 



F = 

According to (21) 



LSA^^^wr- (29) 



OF 



Therefore, taking into account (26) and (29), 

_ 2a(y^) 3 " 
3nb 2 v 11 

Applying the first equation (21), we obtain state equations 

x dF 2a 
A = — 



P dv 3bvp 

Denoting pressure A by p in accordance with established thermodynamical notation, 
we arrive at a more usual form of state equation: 

pv = ~kT. (30) 

Now let f(z) = e' z . Thus, 



m + 2 ./n m + 2 

Hence, | = ( m + 2 ) = ^ and state equation (30) trnasforms into the classical Clapey- 
ron equation: 

pv = nkT . (31) 
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Now assuming that state equations (30) and (31) are identical under any n, we 
can ask the following question. Is it true that frequency function will be of Gibbs 
form, i.e. f(z) = exp(— z)l The answer appears to be negative. Actually (30) 
and (31) are identical if 

a 3n m + 2 n n 

- = - = —, m = 3n-2 = l,4,7,.... 

With account of (27) and (28) these equations take the following form 

/oo 2 p 2 

r m+3 f(j)dr=(m + 2) J r m+1 f(^-)dr. (32) 

Let / be decreasing at infinity faster than any exponential function. Then by part- 
wise integrating we can represent (32) as follows 

for all m+3 = 3n+l = 4, 7, 10, .... If this equality was true for all non-negative m+ 
3, then according to classical momenta theory [18], the expression in the square 
brackets of (33) would be equal to zero. Hence /' + / = and, therefore, / = 
ce~ z , c = const. However, (33) is not valid for the "majority" of integer values 
of m + 3. Hence, it follows that there is an infinite-dimensional space of frequency 
functions dependent on total energy only, which result in classical thermodynamical 
relations for the perfect gas. 

The work has been partially supported by RFBR (grant No. 99-01-0196) and 
INTAS (grant No. 96-0799) foundations. 

References 

[1] G. V. Gibbs. Thermodynamics. Statistical Mechanics. Moscow, Nauka, 1982, 
P. 384. 

[2] V. V. Kozlov. A Constructive Method of Justification of the Theory of Systems 
with Unilateral Constraints. Prikl. Mekh. i Mat., 1988, V. 52, Iss. 6, P. 883-894. 

[3] P. K. Rashevsky. Riemannian Geometry and Tensor Analysis. Moscow, Nauka, 
1967, P. 664. 

[4] N. G. Brain. Asymptotic Methods in Analysis. Moscow, IL, 1961, P. 247. 

[5] A. Ya Hinchin. Mathematical Grounds for Statistical Mechanics. Moscow- 
Leningrad, Gostekhizdat, 1943, P. 128. 



~ m + 3 dr = 



(33) 



12 



[6] V. V. Kozlov. Canonical Gibbs Distribution and Thermodynamics of Mechani- 
cal Systems with a Finite Number of Degree of Freedom. Regular and Chaotic 
Dynamics, 1999, V. 4, No. ,2, P. 44-54. 

[7] E. T. Whittaker, G. Robinson. The Calculus of Observations. Blackil and Son, 
1928. 

[8] A. M. Kachan, Yu. V. Linnik, S. R. Rao. Characterizational Problems of Math- 
ematical Statistics. Moscow, Nauka, 1972, P. 656. 

[9] E. A. Sidorov. Smooth Topological Transitive Systems. Mathematical Notes, 
1968, V. 4, No. 6, P. 751-759. 

[10] J. C. Oxtoby. Note of Transitive Transformations. Proc. Mat. Acad. Sci. U.S., 
1937, V. 23, P. 443-446. 

[11] V. V. Kozlov. Phenomena of Nonintegrability in Hamiltonian Systems. Proc. 
Int. Congr. Math. Berkeley. California. USA, 1987, P. 1161-1170. 

[12] N. G. Moshchevitin. On Existence and Smoothness of an Integral of a Hamil- 
tonian System with a Defined Form. Mathematical Notes, 1991, V. 49, No. 5, 
P. 80-85. 

[13] V. V. Kozlov. Symmetries, Topology and Resonances in Hamiltonian Mechanics. 
Springer- Verlag, 1996, P. 378. 

[14] S. V.BoloUn. Birkhoff Integrable Billiards. Vestnik MGU, Ser. Mat., Mekh., 
1990, No. 2, P. 33-36. 

[15] M. Kac. Probability and Related Topics in Physical Sciences. Intersience Pub- 
lishers, 1957. 

[16] V. F. Lazutkin. A Convex Billiard and Eigenfunctions of the Laplace Operator. 
Leningrad, LGU Publishers, 1981, P. 196. 

[17] N. V. Svanidze. Existence of Invariant Tori for a Three-Dimentional Billiard 
Located in a Neighbourhood of a "Closed Geodesic on the Domain Boundary" . 
UMN, 1978, V. 33, Iss. 4, P. 225-226 

[18] N. I. Ahieser. Classical Problem of Moments. Moscow, Nauka, 1961, P. 310. 



13 



Figure 1: 



14 



